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Abstract 

We have derived the relationship between polarization and intensity fluctuations due 
to point sources. In the case of a Poisson distribution of a population with uniform 
evolution properties and constant polarization degree, polarization fluctuations are 
simply equal to intensity fluctuations times the average polarization degree. Con- 
servative estimates of the polarization degree of the classes of extragalactic sources 
contributing to fluctuations in the frequency ranges covered by the forthcoming 
space missions MAP and Planck Surveyor indicate that extragalactic sources will 
not be a strong limiting factor to measurements of the polarization of the Cosmic 
Microwave Background. 
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1 Introduction 



There are good prospects that the forthcoming space missions designed to pro- 
vide high sensitivity and high resolution maps of the cosmic microwave back- 
ground (CMB) will also measure the CMB polarization fluctuations (Knox 
1998; Bouchet et al. 1999). 
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The current design of instruments for the Planck Surveyor mission (the third 
Medium-sized mission of ESA's Horizon 2000 Scientific Programme) provides 
good sensitivity to polarization at all LFI (Low Frequency Instrument) fre- 
quencies (30, 44, 70, and 100 GHz) as well as at three HFI (High Frequency 
Instrument) frequencies (143, 217 and 545 GHz). The NASA's MIDEX class 
mission MAP has also polarization sensitivity in all channels (30, 40 and 90 
GHz). 

The extraction of the very weak cosmological polarization signal requires both 
a great sensitivity of the instruments and a careful control of foregrounds. An 
analysis of the effect of Galactic polarized emissions (synchrotron and dust) 
on CMB measurements by the Planck and MAP missions was carried out by 
Bouchet et al. (1999). So far, however, the effect of extragalactic sources was 
not considered. On the other hand, significant linear polarization is seen in 
most compact, flat-spectrum radio sources which are the main contributors 
to small scale foreground intensity fluctuations at A > 1 mm (Toffolatti et al. 
1998) and the thermal dust emission from galaxies which dominate the counts 
at sub-mm wavelengths is also expected to be polarized to some extent, as the 
Galactic dust emission is observed to be (Hildebrand 1996). 

In this paper we derive the relationship between intensity and polarization 
fluctuations in the case of a Poisson distribution of point sources, discuss the 
polarization degree of the relevant classes of extragalactic sources, and exploit 
recent evolutionary models to estimate the power spectrum of polarization 
fluctuations produced by them. 



2 Polarization fluctuations from a Poisson distribution of polarized 
point sources 



Following Burn (1966) we define the complex linear polarization of a source as 
P s = nexp(2i%) where n and x are the degree and the angle of polarization, 
respectively. If the polarization angles of different sources within any given 
solid angle element dQ are uncorrelated, the expected value (P s ) is and the 



variance is: 



TV TV 
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Let N = n(S) dS dQ be the number of sources with flux S within dS and 
polarization degree II in a given solid angle element dVL. As far as the central 
limit theorem holds, the expected value of the linear polarization within dVL 



2 



is also 0, with variance 



a p,dn — • (2) 

The fluctuation amplitude of the polarized flux Sp = NSP S among the differ- 
ent cells of the sky subtending a solid angle dfl, due to sources with flux S, is 
therefore obtained integrating over the probability distribution of N, p(N): 

oo 

<4 P <n = ((Sp - (Sp)) 2 ) = (Si) = J dNp(N)NSW = tfS 2 (N(S)). (3) 

o 

In the case of a Poisson distribution of sources the variance is equal to the 
mean. Therefore, integrating the above equation over S and over the solid 
angle we straightforwardly obtain oi p = Hcri, where 07 is the rms intensity 
fluctuation for a Poisson distribution of sources. 

The assumption of an equal polarization degree for all sources is obviously 
unrealistic. However, it follows from the above calculations that, if the polar- 
ization degree is uncorrelated with flux, the result depends only on the mean 
value of n. 

A dependence of the mean value of II on S arises, in particular, in the case of 
contributions from classes of sources with different polarization properties and 
different shapes of the logiV-logS' curves. These different populations must 
be dealt with separately. 



3 Linear polarization properties of the relevant classes of sources 

Studies of the polarization properties of extragalactic sources at mm wave- 
lengths are still scanty. Table 1 summarizes the main results on radio/mm 
sources. For each of the main source classes (in column 1) and for each wave- 
length in column 2, we give the median (column 3), the minimum (column 4) 
and the maximum (column 5) polarization value found in literature. Column 
6 lists the references from which the reported values are derived. 

As shown in Table 1, the three classes of objects (BL Lacs, QSOs and bright 
radio galaxies) emit almost the same fraction of polarized radiation in the radio 
band. The major extragalactic contributors to the non-thermal polarization at 
sub-mm and mm wavelengths are the flat spectrum (spectral index a < 0.5 
if S u oc v~ a ) compact radio sources (mainly BL Lacertae objects and flat 
spectrum QSOs). These objects constitute about 50% of the flux limited (S > 
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Table 1 

Radio and mm polarization measurements for extragalactic sources 







Band 


Pmed 

(%) 


p . 

1 mm 

(%) 


p 

1 max 

(%) 


Reference 


BL Lac 




4.8 GHz 


3.6 


1.5 


7.5 


Aller et al. 1999 






8.44 GHz 


1.5 


0.4 


2.7 


Marcha et al. 1996 






14.5 GHz 


5.0 


1.1 


11.4 


Aller et al. 1999 




0.8 


— 1.1 mm 


7.1 


3.0 


12.8 


Nartallo et al. 1998 




0.8 


— 1.1 mm 


10.8 


3.0 


17.0 


Stevens et al. 1996 


FS QSO 


1.4 


- 90 GHz 


2.5 






Saikia & Salter 1988 


HPQ 


0.8 


— 1.1 mm 


7.4 


4.6 


10.5 


Nartallo et al. 1998 


LPQ 


0.8 


— 1.1 mm 


4.9 


2.1 


7.9 


Nartallo et al. 1998 


FRII 




1.4 GHz 


~4 






Saikia & Salter 1988 






1.4 GHz 


~10 


2.8 


18.4 


Ishwara-Chandra et al. 1998 






5 GHz 


~6 






Saikia & Salter 1988 






5 GHz 


~10 


2.4 


18.2 


Ishwara-Chandra et al. 1998 


FRI 




5 GHz 


~7.5 






Saikia k Salter 1988 



1 Jy at 5 GHz) radio catalogue compiled by Kiihr et al. (1981). Stickel et 
al. (1991) have drawn from this catalogue a complete sample of BL Lacertae 
objects brighter than m = 20 mag. Out of it, we have selected a complete 
sub-sample of 14 objects (RA > 9 h and 5 > 0°), for which radio and/or mm 
polarization measurements are available. 

The polarization data for objects in this sub-sample (hereinafter Stickel north) 
are listed in Table 2. The tabulated percentage polarization degrees were ob- 
tained averaging the measurements from long term monitoring programs. Data 
at 4.8 GHz, 8.1 GHz and 14.5 GHz are from Aller et al. (1985, 1999), those 
at 22 GHz, 31 GHz and 90 GHz are from Rudnick et al. (1985) and those at 
1.1 mm are from Nartallo et al. (1998). 

We have compared the average polarization percentages at the various wave- 
lengths of sources in the Stickel north sample and in the sample by Aller et al. 
(1999), excluding those in common with the Stickel north sample. The agree- 
ment is very good (see Fig. 1). The mean polarization percentages (P m ), the 
68% confidence uncertainties (AP m ) obtained from the Student's t distribu- 
tion (AP m = ^(0.32)^ (N - 1)~ 1/2 , where N is the number of sources in the 
bin and v — N — 1) and the dispersions, ap, for the Stickel north sample are 
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Table 2 

Fractional Polarization Measurements for the Stickel north BL Lac Sam- 
ple at several frequencies (GHz) 



nfl TT1P 


Pf4 8 s ) 


Pf8 1) 


P(14 5) 


Pf22) 




Pf90) 


P(272) 


2 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 




0Q54+658 


7.5 




6.0 










0.367 


1147+245 


2.4 


3.3 


4.6 














2.3 


2.5 


3.0 


5.6 


5.3 


2.0 


9.9 


0.997 


1418+546 


1.9 


2.4 


2.6 






5.3 




0.152 


1538+14Q 


4.9 


5.8 


7.5 










0.605 


1652+398 


1.7 


3.1 


3.0 










0.033 


174g+oqfi 


3.1 


3.3 


3.4 


6.5 




5.5 


5.6 


0.320 


1749+701 


3.6 


5.6 


8.1 










0.770 


1803+784 


3.1 


3.0 


3.9 










0.684 


1807+698 


2.1 


4.1 


2.5 










0.051 


1823+568 


4.3 


5.9 


5.8 










0.664 


2007+777 


3.0 


5.1 


7.1 










0.342 


2200+420 


5.0 


3.6 


5.2 


0.8 


3.2 


4.3 


8.1 


0.069 


2254+074 


6.3 




11.4 










0.190 



given in Table 3. 

Since for most objects in the two samples redshift information is available, we 
have converted the observed into the rest-frame wavelengths and computed 
the average polarization degrees in wavelength bins. The results for the Stickel 
north sample are given in Table 4 and plotted in Figure lb, where the results 
for the Aller sub-sample are also reported for comparison. Table 4 gives the 
adopted wavelength intervals, the average A in each interval, computed as the 
geometric mean of the two wavelength limits, the average polarization fraction 
with its 68% confidence uncertainty from the Student's t distribution and its 
dispersion, and the number of available measurements in each wavelength 
interval. 
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Table 3 

Average Polarization Percentages for Stickel north Sample 



A 




p 


AP 

" J TO 


a P 


N 


(cm) 


(GHz) 


(%) 


(%) 


(%) 




6.25 


4.80 


3.66 


0.49 


1.74 


14 


3.70 


8.10 


3.98 


0.40 


1.29 


12 


2.07 


14.5 


5.29 


0.73 


2.57 


14 


1.36 


22.0 


4.30 


2.83 


3.06 


3 


0.97 


31.0 


4.25 


2.69 


1.48 


2 


0.33 


90.0 


4.28 


1.10 


1.60 


4 


0.11 


272 


7.87 


2.00 


2.16 


3 



Table 4 

Average rest-frame polarization percentage for the Stickel north sample 



AA 


A m 


p 


AP m 


a P 


N 


(cm) 


(cm) 


(%) 


(%) 


(%) 




3.20-6.40 


4.52 


3.65 


0.40 


1.57 


17 


1.60-3.20 


2.26 


4.39 


0.69 


2.43 


14 


0.80-1.60 


1.13 


4.99 


0.69 


2.19 


12 


0.40-0.80 


0.57 


5.45 


0.38 


0.21 


2 


0.20-0.40 


0.28 


5.03 


0.59 


0.64 


3 


0.10-0.20 


0.14 


5.05 


7.83 


4.31 


2 


0.05-0.10 


0.07 


7.75 


5.53 


3.04 


2 
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0.1 1.0 10.0 

Observed Wavelength (cm) 



1 2 




o r i , , , , i 

0.1 1.0 10.0 

Rest Frame Wavelength (cm) 

Fig. 1. Average fractional polarization as a function of wavelength for two sam- 
ples of BL Lac objects: red filled circles connected by the red dashed line stand for 
the Stickel north sample, while the blue filled stars connected by the blue dotted 
line stand for the the 28 sources in the Aller et al. (1999) sample that are not in- 
cluded in the Stickel north sample. The error bars correspond to the 68% confidence 
uncertainties on the estimated mean polarization. 
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Measurements of polarized thermal emission from dust are only available for 
interstellar clouds in our own Galaxy. The distribution of observed polarization 
degrees of dense clouds at 100 /im shows a peak at ~ 2% (Hildebrand 1996). 
The polarization degree of emission of silicate grains in clouds opaque to visible 
light but optically thin in the far-IR is nearly independent of wavelength if 
A ^> a, a being the grain size (Hildebrand 1988). This condition is very likely 
to be met at the long wavelengths of interest here. 

Polarization maps of the Orion molecular cloud at 100 /im, 350 /im, 450 /im, 
1.3 mm, and 3.3 mm (Schleuning 1998; Rao et al. 1998) look very similar ex- 
cept at clumps of higher optical depth, where the polarization increases with 
wavelength. In fact, the maximum polarization decreases rapidly with increas- 
ing optical depth (Hildebrand 1996); thus, the mm/sub-mm polarization may 
be higher than at 100 /im. On the other hand, the overall polarization of the 
light from a galaxy is the average of contributions from regions with different 
polarizing efficiencies and different orientations of the magnetic field with re- 
spect to the plane of the sky; all this works to decrease the polarization level 
in comparison with the mean of individual clouds. 

We could not find any published polarization measurement of dust emission in 
external galaxies not hosting strong nuclear activity. However, the first results 
of the SCUBA Polarimeter include imaging polarimetry of the starburst galaxy 
M 82 at 850 /im with a 15" beam size; a polarization degree of about 2% was 
measured (J. Greaves, private communication). An early image, available at 
the polarimeter Web page ( [http: / / www .jach.hawaii.edu/ J UMT/ scuba/ scupol/ 



general/m82_polmap.gif), shows very ordered polarization vectors over scales 
of a few hundred parsecs. This may not occur in general. We may expect 
that, in other galaxies, the polarization vectors from giant molecular clouds 
are randomly ordered and tend to cancel out. 

Polarimetric measurements at 170 /im with ISOPHOT have been carried out 
for two galaxies (U. Klaas, private communication): NGC 1808 (P.I.: E. Kriigel) 
and NGC 6946 (P.I.: G. Bower). The results, however, are not available yet. 



4 Power spectrum of polarization fluctuations due to extragalactic 
sources 



In Figures 2 and 3 the power spectrum of foreground polarization fluctuations 
for all the relevant Planck and MAP channels is compared with the power 
spectrum of CMB anisotropics and of CMB polarized components. As for the 
latter, we have plotted the power spectra of the combinations of Stokes pa- 
rameters defined by Seljak [1997; his eqs. (24) and (25)] and called E and B. 
The estimate of E'-mode polarization fluctuations refer to a standard CDM 
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30 GHz 



40-44 GHz 



70 GHz 








000 



100 
Multipole (/) 

Fig. 2. Power spectra of CMB brightness temperature (thick solid line) and polar- 
ization fluctuations (.E-mode: long dashes; B-mode: three dots/dash) compared with 
foreground polarization fluctuations at Planck/LFI frequencies. Following Tegmark 
& Efstathiou (1996) we plot the quantity 5T e (u) = [£(2e+l)C t (v)/4:ir\ 1 ^. The thick 
dot-dashed and dashed lines correspond to synchrotron and interstellar dust emis- 
sions, respectively. The thin solid and dashed lines show our estimated contributions 
from extragalactic radio and far-IR sources, respectively, based on the model by Tof- 
folatti et al. (1998), as updated by De Zotti & Toffolatti (1999). The thin dotted 
lines represent the expected noise spectrum per resolution element, averaged over 
the sky, of Planck/LFI. The heavier dots show the average MAP noise spectrum 
per resolution element at 30, 40 and 90 GHz. 
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10 100 1000 
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Fig. 3. Power spectra of polarized components at Planck/HFI frequencies with 
sensitivity to polarization. The lines have the same meaning as in Fig. 2. We show 
estimates of polarization fluctuations due to dusty galaxies (thin dashed lines) based 
both on model E by Guiderdoni et al. (1998; upper line) and on the model by 
Toffolatti et al. (1998). 

model (scale-invariant scalar fluctuations in a O = 1, A = universe with 
H = 50kms _1 Mpc -1 and a baryon density = 0.05). The quantity B van- 
ishes for polarization induced by primordial scalar perturbations, and therefore 
provides a unique signature of tensor perturbations (Seljak 1997). The I?-mode 
power spectrum shown in Figs. 2 and 3 refer to a tilted CDM model with power 
law indices n s = 0.9 and n t = n s — 1 = —0.1 for scalar and tensor pertur- 
bations, respectively; the other cosmological parameters keep the same values 
adopted for the standard CDM model. Calculations of CMB power spectra 
have been carried out using the CMBFAST package by Seljak & Zaldarriaga 
(1996). 

The thick dashed lines show the i?-mode of dust polarized power spectrum 
derived by Prunet et al. (1998), scaled to the central frequencies of Planck 
polarized channels using the dust emission spectrum adopted by these Authors 
(emission oc v 2 B v (17.5 K), where B V (T) is the Planck function at frequency v 
and temperature T). 

The .B-mode dust power spectrum turns out to be close to the E-mode one 
[cf. eqs (11) and (12) of Prunet et al. (1998)]. In fact, Seljak (1997) argued 
that most foregrounds should contribute on average the same amount to both 
modes. Therefore, for all foregrounds we consider a single polarized power 
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spectrum, assumed to be representative of both modes. 

Following Bouchet et al. (1999), we assume that, at all frequencies relevant 
for Planck and MAP, the synchrotron polarized emission is perfectly cor- 
related with the total synchrotron emission (for which a power spectrum 
Ct = 4.5£~ 3 (/iK) 2 at 100 GHz was adopted) and the polarization degree 
is 44%. An antenna temperature spectral index of 3 (TA, sy n oc z/~ 3 ) has been 
used to extrapolate the power spectrum (in terms of brightness temperature) 
to the other Planck/MAP frequencies (see De Zotti et al. 1999 for references). 

The polarized angular power spectra of extragalactic radiosources (thin solid 
lines) is estimated exploiting the model of Toffolatti et al. (1998), as updated 
by De Zotti & Toffolatti (1999). We have assumed that the mean polarization 
degree of BL Lacs in the Stickel north sample applies to all radio galaxies 
contributing to fluctuations in Planck's channels; based on the results shown 
in Fig. 1, we have adopted a polarization degree of 5% for v < 143 GHz, of 
6% at 217 GHz and of 10% at 545 GHz. 

As for dusty galaxies (thin dashed lines), we have adopted a polarization degree 
of 2% at all frequencies and the power spectra of temperature fluctuations 
derived by Toffolatti et al. (1998) and, for HFI channels, also by Guiderdoni 
et al. (1998; their model E). 

A flux cut at 1 Jy was adopted for all channels (i.e. sources brighter than 1 
Jy were removed); at 30, 100 and 217 GHz we show also the angular power 
spectrum derived adopting a flux cut of 100 mJy for radiosources (lower thin 
solid lines). 

In Table 5 we report the values of Cg for temperature fluctuations due to a 
Poisson distribution of extragalactic radio sources and dusty galaxies for the 
same cases shown in Figs. 2 and 3. The corresponding values for polarization 
fluctuations follow immediately multiplying by n 2 . Note that a Poisson dis- 
tribution generates a simple white noise spectrum with the same power on all 
multipoles (Tegmark & Efstathiou 1996). The values of Ci are given in terms 
of brightness temperature fluctuations and expressed in fiK 2 . 

As expected, polarization fluctuations due to extragalactic sources are partic- 
ularly relevant at small angular scales. For multipoles i ^ 300, they in fact 
dominate foreground contributions at v ;$ 100 GHz. On these scales, at the 
lowest Planck frequency (30 GHz) their amplitude is, according to our esti- 
mate, close to that of CMB polarization fluctuations induced by scalar per- 
turbations. In the "cosmological window" (70 ;$ v £j 200 GHz), however, ex- 
tragalactic sources are not seriously detrimental to measurements of CMB 
polarization fluctuations. 
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Table 5 

Values of Cg (fiK 2 ) 1 for temperature fluctuations for the cases shown in Figs. 2 and 
3 



v (GHz) 



30 



44 



70 



100 



143 



217 



545 



Si = Uy 

Radiosources 5.0(-2) 1.5(-2) 3.2(-3) 0.9(-3) 2.5(-4) 9.5(-5) 5.6(-3) 
(Toffolatti) 

Far-IR sources 1.5(-4) 4.0(-5) 7.6(-6) 3.1(-6) 2.8(-6) 1.0(-5) 3.0(-l) 
(Toffolatti) 

Far-IR sources 5.4(-5) 2.5(-4) 8.8(-l) 

(Guiderdoni) 

Si = 100 mJy 

Radiosources 1.0(-2) 1.0(-4) L2(-5) 
(Toffolatti) 

Far-IR sources 5.0(-5) 1.5(-6) 6.0(-6) 
(Toffolatti) 

Far-IR sources 2.0(— 4) 
(Guiderdoni) 

1 In parenthesis are the powers of 10 (i.e. 5.0(— 2) = 5.0 10 -2 ) 

The thin dotted lines in Figs. 2 and 3 show the expected power spectra of 
instrumental noise, for polarization measurements, averaged over the sky, for 
Planck's LFI and HFI, respectively. Following Tegmark & Efstathiou (1996) 
we describe the noise power spectrum as C^ no i SC = <r 2 FWHM 2 exp(£ 2 a 2 ) where 



FWHM in expressed in radians, <7& = FWHM/2y2 ln(2) and a is the rms noise 
for a square pixel with side FWHM. 

For LFI we have adopted the sensitivities for brightness temperature measure- 
ments given by Mandolesi et al. (1998) multiplied by a factor of 2 (Mandolesi, 
private communication). Sensitivities of HFI channels for polarization mea- 
surements are given by Puget et al. (1998). There is a slight difference in the 
mission duration adopted by the two groups to derive their mean sensitivity 
estimates: 12 months for LFI, 14 months for HFI. 

The heavier dots in Fig. 2 show the expected mean instrumental noise per res- 
olution element for MAP's polarization measurements at 30, 40 and 90 GHz 
obtained from sensitivities for brightness temperature measurements (see the 
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MAP Web page) multiplied by a factor \pl (G. Hinshaw, private communica- 
tion) . 

As shown by Figs. 2 and 3, the major hurdle in extracting the CMB polar- 
ization signal is instrumental noise. However, a simple argument shows that, 
at least for a limited band-power range, sufficient sensitivity can be reached. 
We have investigated, in particular, the potential of LFI in this respect. The 
expected sensitivities to polarization per resolution element, averaged over the 
sky, for a 12 months mission, are 6, 10, 14 and 17/iK at 30, 44, 70, and 100 
GHz, respectively; the angular resolutions (FWHM) are 33', 23', 14', and 10', 
respectively (Mandolesi et al. 1998). Simulations done by C. Burigana indicate 
sensitivities 7 times better over areas of about 25 square degrees around each 
of the ecliptic poles. Within these areas, by rebinning the maps at 44, 70 and 
100 GHz (we leave aside the 30 GHz map which is the most contaminated by 
polarized foregrounds) to a ~ 20' resolution and combining them, a sensitivity 
to polarization of about 1 /iK can be achieved, allowing to image the CMB 
polarization induced by scalar perturbations predicted by the standard CDM 
model. Of course, a lower sensitivity is enough to determine the first moments 
of the distribution of polarization fluctuations. 



5 Conclusions 

We have shown that the polarization fluctuations due to a Poisson distribu- 
tion of point sources with uniform evolutionary properties and constant po- 
larization degree are simply equal to intensity fluctuations times the average 
polarization degree. 

The information on the polarization degree of the classes of extragalactic 
sources expected to dominate in the frequency ranges relevant for the MAP 
and Planck missions is scanty. We have taken a conservative approach, as- 
suming that all radio sources are as polarized as BL Lac objects and that the 
polarization degree of dusty galaxies is similar to that of dense clouds in our 
own Galaxy and of M82, a galaxy showing a remarkably ordered magnetic 
field. 

We find that, on small scales (multipoles £^,300), polarization fluctua- 
tions due to radio sources may indeed dominate foreground contributions at 
v< 100 GHz. However, in the "cosmological window" (70 < v £ 200 GHz), 
extragalactic sources are not a threat for measurements of CMB polarization 
fluctuations. 

We have also argued that Planck/LFI can reach, in regions around the Galac- 
tic polar caps, polarization sensitivities of ~ 1 /iK, allowing to map CMB 
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polarization fluctuations on scales ~ 20'. A detailed analysis of Planck and 
MAP capabilities for CMB polarization measurements has been carried out 
by Bouchet et al. (1999). 
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